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1. INTRODUCTION 


The primary goal of this paper is to illustrate, by a simple problem, the 
necessity of conducting a careful analysis of numerical schemes (that were de- 
veloped for open-loop simulation) when these schemes are to be applied to an 
optimization based control design problem. It is our feeling that many dis- 
tributed parameter control systems (viscoelastic structures, fluid flow con- 
trol, etc.) are such that "standard finite element/finite difference" schemes 
might lead to numerical difficulties in certain control design methods, if 
care is not taken to ensure that these approximation schemes have convergence 
properties essential for the intended application. 

We shall concentrate on a non-convergence result for an optimal control 
problem governed by a delay differential equation. Although we feel that the 
technical details required to prove non-convergence are interesting, we hope 
that this proof is not viewed as the major contribution of the paper. Indeed, 
we hope that the reader is motivated to think about similar problems for more 
complex distributed parameter systems. 

During the past fifteen years considerable attention has been devoted to 
the construction of finite dimensional approximations of distributed parameter 
systems. Much of this work is based on algorithms first developed primarily 
for simulation. However, it is not clear that finite dimensional models 
developed for (open-loop) simulations will also be suitable for certain opti- 
mization based control design techniques. Moreover, there may be several 
reasons that a numerical scheme developed specifically for simulation does not 
perform well when applied to a control design problem. 

In this note, we concentrate on the use of the so-called spline scheme 
developed by Banks and Kappel [2] as an approximation method for calculating 



optimal feedback gains for control systems governed by functional differential 
equations. This problem is simple enough to be addressed in a short paper and 
yet it still illustrates how a specific numerical scheme when combined with a 
particular control design approach can lead to numerical difficulties. 

In order to motivate the paper and to describe the main technical con- 
tribution of the paper, we first review some known results. 

Let H and U be real Hilbert spaces and S(t): H ■*> H denote a Cg- 
semigroup of bounded linear operators with generator A. We assume that 
B: U H, Q: H -*■ H and R: U ->• U are bounded linear operators with Q 

and R self-adjoint, non-negative and R satisfies R ^ ml > 0. The (LQ) 

optimal control problem is to minimize 

T 

(1.1) J(u) = / [ <Qz(s) , z(s)> + <Ru(s) ,u(s)>]ds 

0 

where z(t) is defined by 

t 

(1.2) z ( t ) = S(t)z + / S(t-s)Bu(s)ds 

0 

for 0 _< t < T and given zq € H. If T = +°°, then one has the linear 
quadratic regulator problem. 

Assume now that there exists a sequence of Cg-semigroups S^(t) on H 
and positive constants M, £5 such that 

(1.3) IIS N (t)ll < Me 6t , t > 0, N > 1 


and 
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(1.4) S N (t) S(t) strongly as N + +» 

and the convergence in (1.4) is uniform in t on bounded intervals. Denote 

by the generator of S^(t) and also assume the existence of bounded 

N N m 

linear operators B: U H, Q : H + H with Q w self-adjoint and 

_> 0 satisfying 

(1.5) »B N - B« + 0, I Q N - Q# -► 0. 

Note the assumption of uniform convergence (1.5) is stronger than required by 
Gibson (see page 113 in [4]). It is well known that the optimal control (if 
it exists) is given by state feedback 

(1.6) u*(t) = -R -1 B*n(t)z*(t) 

where the bounded self-adjoint operator II ( • ) satisfies the Riccati 

(integral) equation 

T 

(1.7) n(t)z » J S*(s - t) [Q - n(s)BR _1 B*n(s) ]S(s-t)ds. 

t 

N 

Let II (t) be the solution to the "approximating" Riccati equation 

(1.8) n N (t)z = / S N *(s-t)tQ N - n N (s)B N R~ 1 B N *n N (s)]S N (s-t)ds 

t 

and observe that (1.8) would determine the sub-optimal gains if one used the 
approximating system (A^,B^) with weights Q N . 
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The following theorem is a direct consequent of Gibson's (more general) 
results (Theorem 6.1 and Theorem 6.2 in [4]). 

Theorem 1.1: If conditions (1.1) ~ (1.5) hold, then for 0 <, t <_ T, 

n N (t) converges weakly to n ( t ) and the convergence is uniform on 

[0,T]. If in addition 

(1.9) S N *(t) ->• S*(t) strongly, 

then II N (t ) converges strongly to n ( t ) and the convergence is uniform 

for t e [0,T] . 

Corollary 1.2: Let_ K(t) = R -1 B*n(t) and K N (t) = R -1 B N *n N (t) 

denote the feedback gain operators. Assume that U = IE? (i.e., is finite 
dimensional). If II^(t) converges strongly to II ( t ) , then as N -*■ +«, 

(1.10) II K N (t) - K(t)H > 0. 

The main point to be emphasized is that if the control space is finite 
dimensional, then uniform convergence is assured provided the numerical scheme 
is stable and consistent and (1.9) holds. If one is concerned only with simu- 
lation, then stability and consistency is sufficient for most numerical prob- 
lems. Moreover, it can be shown that many standard numerical schemes de- 
veloped for simulation of self-adjoint partial differential equations do 
satisfy (1.9). Therefore, it is not surprising that until Gibson "needed” 
(1.9) to establish the uniform convergence of optimal gain operators, the 
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question of whether a numerical scheme satisfied (1.9) received little atten- 
tion. Indeed, even after Gibson published his result it was still not obvious 
that condition (1.9) was anything more than a technical assumption needed in 
Gibson's proof. 

In [3] Banks, Ito, and Rosen applied a convergent spline based scheme to 
an optimal control problem governed by a delay differential equation. The 
numerical results in [3] seemed to show that K^(t) did not converge uniform- 
ly to K(t) and these numerical results have often been used as evidence that 
(1.9) did not hold for this particular scheme. Moreover, several new schemes 
have since been generated specifically to ensure (1.3) - (1.4) and (1.9) are 
valid. Still, it was not known if condition (1.9) held for the spline scheme 
used in [3]. We shall provide a proof that (1.9) fails for this scheme. We 
also show that this spline scheme is stable and consistent to A on a dense 
subset of D(A ). 

2. SPLINE APPROXIMATIONS OF HEREDITARY SYSTEMS 

Consider the delay differential equation 

0 

(2.1) x(t) = ApX(t) + AjX(t - r ) + / A(s)x(t + s)ds 
with initial data 

(2.2) x(0) = n ; x(s) = <|>(s), -r < s < 0, 

where x(t) I? and the elements of A(») are square integrable on 
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[ — r » 0 ] • It is well known that (e.g., [1]) for n € If 1 and 

<j> e L^(-r,0; I?) there exists a unique solution of (2.1) - (2.2) 

x: [-r,+°°) I? such that x€W^’^(0,T; I?) for all T > 0. If one 

defines the solution map S(t), t > 0 on the product space 

Z = IP x L^-rjO;!? 1 ) by 

(2.3) s(t)(n ,<(>(• )) = (x(t), x(t+»)) 

where x is the solution to (2.1) - (2.2), then { S( t ) } a stron S^y 
continuous semigroup (i.e., C 0 -semi g rou P) on z * The infinitesimal 
generator A is the operator defined on the domain 

(2.4) D(A) = {(p , 4 > (- )) ez|4,(.)ew 1 * 2 (-r,0; If 1 ), 4 . (0) = n} 
by 

0 

(2.5) A(n ,$(•)) = (A q ti + A 1 <j»(-r) + / A(s)<{> (s)ds , <(>(•)). 

-r 

The adjoint operator A* generates the adjoint semigroup S*(t) and it is 
easy to show that (see [4,7]) 

(2.6) D( A*) = {(? ,i|>) £ Z|i|) € W 1,2 (-r,0; If 1 ), ^(-r) = A^} 
and for (g ,i|/) € D(A*) 

A *(?,*) = 0K0) + aJc, [A T (.)5 - *(•)])• 


(2.7) 
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i 


! 


I 


i 


As in [ 2 ] , 

Let 

[-r,0] by 


(B N 1 N 
{B i } i«0 


we define the linear spline based approximation for S(t)« 
denote the usual linear B-splines defined on the interval 



j 

N 

(S - T 1+1 ), 

se [x N x N ] 
6 lT i+i ,T i J 


| 

r 

(2.8) 

Bj(s) = | 

‘ N 

r N x 

(T i-i - s) - 

- r N N 

selt l- T i-1 



0 

y 

otherwise, 


where = -ir/N, i = 0,1, = -r and = 0. For each 

N = l 9 2 9 ^^^ let Z N denote the linear subspace of Z defined by 

N 

(2.9) Z N = (ze z|z = l \(b^(0),b^(.)), a^iP} 

k=0 

and let denote the orthogonal projection of Z onto Z^* This subspace 

can be identified with ]jp(N+l) by the prolongation i^s ]gi(N+l) + £ 
defined by 

XT N M 

(2.10) ia = (a , £ iB N (.)) 

U k=0 

where a = (aL aT,» • • .a^^eiP^+l) . The space is norraed with the 

0 1 N 

induced inner product 


(2.11) <a,b> N - a T Q N b, 

where a, b € lP (N+1) and Q N is defined by 
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(2.12) Q N = ° 


—I —I 

r 6 3 • 

+ N * • 


The adjoint operator [i N ]*: Z -► is given by 


(2.13) 


[i N ]*(n ,<(.(•)) = tQ N ] 1 


n + <f>, 


u N 

where <J> ^ = / 4>(s)B;;(s)ds. Moreover, it is easy to show that 


(2.14) 


[i N ]*i N - I the identity on lf l(N+1) 


i N [iV - P N 


and for z, w£ Z 


(2.15) 


<z,w> = <[i N ] z, [i N ] w > N< 


In order to construct the standard Galerkin approximation of A, we note 
that Z N C D( A) and define A N by A N = P N AP N . Observe that A N (and 
hence [A N ]*) is continuous and although P N Z C D(A), P N does not map all 
of Z into D( A*) . It is shown in [2] that 


a n = i n [ Q ■rVu 11 ]* 


(2.16) 


where 


(2.17) 


and 


1 

• 

• 

• 

+ 1 

’l - 1 . 

I 0 

0 


• • 

• 

• 

0 

2 

• 0 

-I 



0 * I 

-I 


Aq = A q + / A(s)Bg(s)ds 

-r 

N 0 N 

A k - / A(s)B”(s)ds, 1 <. k < N-l 

A^ = A x + / A(s)B^(s)ds. 


If S«(t) denotes the CQ-semigroup generated by A^, then it is shown in [2] 
that for each z€Z, 


(2.18) 


IS N (t)z - S(t)z! 0, 


where the convergence is uniform in t on bounded intervals. The convergence 
(2.18) was established by proving that the Galerkin approximations A N 
satisfy the Trotter-Kato Theorem [6] and hence provided a stable and consis- 
tent approximation scheme for A. We shall prove that the above convergence 
statement does not hold for the sequence of adjoint semigroups [S^]*(t). 
Moreover, this convergence fails even though the adjoint operators [A^]* are 
stable and consistent to A* on a dense subset of Z (i.e., there is a dense 
set C C D(A*) c Z such that [A N ]*z A*z for all z € C) . 
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3. CONVERGENCE OF THE ADJOINT GENERATORS 

We shall follow the approach given in [5] and consider the n(N+l) x n 
matrix 

N r „N, — 1 f e\ n t nT 

Y + - [Q ] (0 0,1) . 

Define the operators 6 + : I? by 

6+x = (t x (0), + x (.)) 

where 

4> x ( * ) - [Bq(-)I, bJ(.)I, • ••, B N N (.)ihJx. 


It follows that for all (ri ,<J> ) « z € Z and x € I? 


rr.N,* N ,T N r Ni* N r N ,1^.1*, 
[[i ] 6 + x] Q [i ] P z - [y + x] Q [i ] 


* 

z 


T N N N * 

- x^B^-Ol^.-.BjJC-r)!]!^] z, 


and (2.14) - (2.15) implies 

(3.1) <6^x,z> z = x T [<}> N (-r) ] , 

where P^z = P N (n,<f>) = (<|> N (0) ,<(> N (* )) € Z N . Furthermore, if denotes 

the smallest eigenvalue of 0 N , then (r/6N) < A^. and it follows from 
(2.11) that 


_< (6N/r) 1/2 


(3.2) 


for all N 
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We also need the following representation. 

w * M 

Lemma 3.1: The operators [A ] : Z -► Z are given by 

[A N ]*(C,ilO = PV(0) + Aj/(0), [A T (« )4» N (0) -$ N (*>]> 

+ 6^/(0) - /(-r)), 

where P N (5,'I») = 0|» N (0) ,A* ) ) . 

Proof : Assume that z = (n ,<(> ) and w = (5,i|») belong to Z and 

let p N z = (4> N (0) ,4> N ( • ) ) and p N w = (lAo) ,<A*)) denote the orthogonal 

projections of z and w, respectively. The identity (3.1) implies that 

<(♦"(0) + aJAo), [A T (.)Ao) - A*)]), P N z> z 
+ <6^(A 1 t|, N (0) - A-r)),z> z 

= t/(0)] T [A0)] + [A0)] T A 0 [<fr N (O)] -/ 4 N (s),<|> N (s)>ds 

-r 

+ [A0)] T J° A(sH N (s)ds + [♦■(0» T A 1 [A-r)l 
-r 

N T N 

- [* (-r)] 4 [* C-r>] . 

Integrating by parts, the boundary terms cancel with the first and last term 
in the sum. Therefore, it follows that if [A N ]* is defined as above, then 
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, r .N,* v Xtj N a N . . 4 N . 

<[A ] w,z> = <P w,AP z> = <w, A z> 

and this completes the proof. 

It should be noted that [A N ]*w = P N A*P N w if and only if 
A^i|) N ( 0) - ^ N (-r> = 0, i.e., if and only if P N w€D(A*). Also, if 

D = D(A)Pl D(A*) 

and 

C = {(5 , 1 ))) e D|il» € C 2 (-r,0; I? 1 )} 

then D and C are dense in Z. Moreover, we have the following convergence 
result. 

Lemma 3.2: C and D are defined as above, then 

(a) [A N ]*w -► A*w for all w€C 

and 

(b) for all X £ E, (XI - A )D is not dense in Z. 

Proof : Let w £ C C D( A) O D(A*) . Note that w = (ij;(0) ,i|>(» )) where 

<K~r) " A^ij)(0) and P N w = (t|> N (0) ,i|> N (* )) £ D(A) . It follows from (2.7) and 

Lemma 3.1 that 

II [A N ]*w - A*wll _< IIP N A*w - A*wll 

+ i(* N (0) + aJ/(o),[a t (.)iAo) - * M (.)1) - (*<0) + aJ*(0),[aJ(.)*(0) - $(.)]) 
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N T N N 

+ I«"[aJ+"( 0) - * (-r)]l 


Since w = (*(0) , *(*))£CCD(A)Od(A ) and *(• ) € C 2 (-r ,0; I?) , it 
follows from standard estimates on interpolating splines (see equations (4.1) 
- (4.3) in [2]) that 

k N (0) - *(0)1 < 0(1/N 2 ) 


|/(s) - *(s)| < 0(1/N), -r < s 0 

and 

ll/(* ) “ *(•)!! < 0(1/N). 


The first term F^ -► 0 since w£D(A*) and p^z -»• z for all z€ Z. The 
second term is estimated by 


< 1/(0) + Aj* N (0) - *(0) - aJ*( 0)| + IIA T (.)(* N (0) - *(0))ll 


+ »/(• ) “ *(• )» 

< (1 + |A 0 | + #A(.)H)k N (0) - *(0)| + 0(1/N). 

N 

Therefore, -*■ 0 as N ». Applying (3.2), the last term is estimated 
by 

f” < (6N/r) 1/2 | A^* (0) - * N (-r)| 
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< (6N/r) 1 ( | A 1 1 |/(0) -ip(0)| + | (— r) -/(-r)|), 

and hence + 0 which establishes part (a) . 

Turning to part (b) , let w = (i|^(0) ,i|j) € D, x £ I? and X € 
straightforward calculation yields that 

<(g,e X ’x), [XI - A*] (*(<» »*("))> 

= /(0)[(X - 1)1 - A Q ]g + * T (0)[I - e^'Aj - / e Xs A T (s)ds]x. 

-r 

If (X - 1) € o(A q ) and e Xr €o(A 1 +J e (X+r)s A T (s)ds) , then 
exist g * 0 and x * 0 such that 

<(£ ,e X *x) , [XI - A ] (xp (0) ))> = 0. 

If (X — 1) t a(Ag), let x * 0 and define 5 by 

5 = -[(X - 1)1 - A q ] _1 (I - e _X r Aj - / e Xs A T (s)ds)x, 

-r 

or if e Xr i a(A 1 + /° e (X+r)s A T (s)ds) , let g * 0 and define x 

-r 

x = -[I - e" Xr A 1 - J° e Xs A T (s)ds] -1 [(X - 1)1 - A Q ]g. 

-r 

In any case, there always exists an element z = (£,\|>(»)) = (C»e X x 

<z , (XI - A*)w> = 0 


E. A 


there 


by 


* 0 , 


such that 
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for all w = (ip(0) , i|)(* )) € D. Consequently, (XI - A*)D is not dense in 

Z. 

Since II [S^] ( t ) II = IIS^(t)ll < e u ^ for some w independent on N 
(see [2]), the existence of a set E C D(A*) and a X € H such that 
(XI - A )E is dense in Z and [A w ] w A w for w € E would imply 
strong convergence of the semigroups [S N ]*(t) to S*(t) ([6], III, Th. 

4.5). Although [A N ]*w A*w on the dense set C, (XI - A*)C is not 

dense in Z. We shall establish that there does not exist a set E c D(A*) 
with the above properties. In fact, we shall show that [S^]*(t) does not 
converge strongly to S (t)* 


4. CALCULATION OF P**( [A**]*) -1 P** 

In this section, we present several technical lemmas that will be needed 
in Section 5. The proof of Lemma 4.1 is straightforward and hence omitted. 
Lemma 4.3 can be found in [4,7]. 

Lemma 4.1; The operator [A^]* can be decomposed as 

(4.1) [A N ]* = i N [Q N ] -1 [H N ] T [i N ]*. 

Lemma 4.2: Assume that is invertible. If z = (p ,<)»(•))€ Z, 

then pN( [A^]*)“lpNz exists and 


p N( [A N]*)- 1 pN z = i N a N > 


(4.2) 
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where is the solution of 


(4.3) 


[H N ] T a N - 


n + <f> 

+ 1 

• XT 


Proof : Consider the equation A N * = P N z for w € Z N . By Lemma 4.1 and 

from (2.14) 

. N, Nx— 1„NT .N* ,N,N* 

i (Q ) H iw=iiz. 

N* 

Multiplying with i and from (2.13) and (2.14) 


U NT . N* 

H i w = 


n + <|> 

* N 

*1 


* N 


M M M* M \T 

The lemma now follows from the fact that w = P w = i i w = i a . 


Lemma 4. 


3s If_ A » Aq + Aj + / A(s)ds, then Oep(A*) if and 


-r 


only if A is invertible and (£,*/>(•)) - (A ) ^ (r) ,(}>(•)) 6 D(A ) is given 
by 


(4.4) 


5 = [A T ] 1 (q + J <|>(s)ds) 
-r 


0 


(4.5) 


s s 

<Ks) * (a^ + J A T (e)de )5 - J <t»(e)de 

-r -r 


Consider equation (4.3). It follows from (2.17) that (4.3) is equivalent 


to the system 
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(4.6) 


1 r, N N A r A N,T N a N ...... 

2 U k+1 " a k-l } + [A k ] a 0 = *k 1 £ k < N_1 


(4.7) 


1 , N . N 


1 ( ^ 4 . a « \ . r.N.T N , N 

2 U N + a N-l } + [ V a 0 = ♦n 


(4.8) 


1 /■ a N . N. . r.N.T N _ N 

2 ^ 0 + P + IV a 0 ~ n + +0* 


A straightforward induction argument yields 


(4.9) 


N N o V /TA N ,T N N . . „ 

a 2k a 0 ” 2 1 ^ 1 ^ A 2i-1 a 0 “ *21-1^* if 2k £ N 


and 


(4.10) 


N N _ f /fA N ,T N N > IC „ , , , „ 

a 2k+l a l 2 ^ ^ A 2i a 0 ^2i)’ if 2k + 1 £ N » 


Thus, it follows from (4.7) that 


(4.11) 


1 , N N w r r.N.T N _ ? N 

" 7 (a 0 " a l ) + J, [ V a 0 " J, V 

k=l k=l 


Moreover, (4.8) and (4.11) imply that 


where 


N H t » N 0 

l [A, ] a^ = n+ l <l>. = n+ J <(>(s)ds 
c=0 K U j=0 J -r 


l [aV = [A. + A. + J° A(s)ds] X = A T . 
k=0 -r 


* f 

If one assumes that 0€p(A ), then by Lemma 4.3 it follows that A is 
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invertible and 


(4.12) 


q = A T (n + / <(> ( s )ds ) i 


-r 


N 


Observe that (4.12) implies that a Q is independent of N and by (4.4) it 
follows that a Q = 5 where (C »4» ( * ) ) “ [A*] -1 (n ,<{>(•)) . Equation (4.8) 
yields the identity 


(4.13) 


N _ N 0 ? , r .N,T N , N> 

a l - ~ a 0 + 2 I ( tA k ] a Q <)> k ) 


k=l 


and hence it now follows from (4,10) that 


N 

1 2k+l 


(4.14) 


N ^ _ ? /r* N i T N .N. , r . r .N ,T N N . 

= ~ a 0 + 2 l, ([A k ] a 0 ‘ " 2 / ( A 2i a 0 ' <f> 2i ) 

k=l 1=1 

N „ ? /r.N.T N .N. . „ r . r .N .T N N . 

= -a Q + 2 l ([A ] a Q - 4 ) + 2 \ ( t A 2i _ 1 3 a Q - ♦jj.j). 

U i=2k+l 2 V 2 i=l 


Note that for 1 < k < N - 1 


N 


Ujl 2 = 1/ 4>(s)B k (s)ds| 2 < (fe / |<f>(s) 1 2 ds 


-r 


N 

r k+l 


UqI 2 i ( 3 n } / |<f>(s)| 2 ds 

N 
T 1 


l4>Sl 2 < <3n> / 


N 

r N-l 


-r 


and similarly 
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Therefore, (4.9) and (4.10) imply 

1*2*1 i ^ (1 + 2 l A ll ) + (|IAII L l a ()l + ■♦ l L ) 

and 

I a 2k+l I 1 l a l 1 + 2 1 A 1 1 I4l + 2/773 (IA. |.*| + 114.11^), 
respectively. The identity (4.12) yields the estimate 

I Sq| < !a- t | (1 + /r)l(n ,<()(• ))n z 
and (4.13) leads to the bound 

|aj| * |-(1 + 2[A^] T )a£ + 2(n + ♦ £) | 

< | Sq I + 2 In | + 2/r73N (iai l Uq| + II<J>I L )• 

Combining these estimates one has that for 0 _< k N 

|a£| < Mil (ri , 4> ( - ) ) II z 

where M > 0 is independent of N. An application of Lemma 4.2 yields the 
estimate 

IIP N [A N *]'' 1 P N ZII z = lli N a N H < (/l+r)M II zll z 
for all z£ Z. We summarize these results in the following theorem. 
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Theorem 4.4: If 0€p(A), then 


pWVc.+o) - A" 


where a N - col( [a”] T , [a"] T ,- • • , [ajjl T ) is given by 


= A T (n + / <|>(s)ds), 

-r 


and for 0 < 2k < N 


N N 
a 2k ' a 0 “ 


. v /r»N ,TN .N . 
2 i-1 ^ A 2i“l 3 0 ” *2i-l* 


while for 1 < 2k+l < N 


N N L . ? , f .N,T N ,N> , r n .N ,T H N 

2k+l ‘ - a 0 + 2 <( V S 0 - V + 2 <IA 2i-l’ *0 ' *21- 


Moreover, II P N [ A N * ] ~ 1 P N II ^ ( z ) is uniformly bounded in N. 


Let p“ denote the "averaging" orthogonal projection on Z 

AVE 

by 

N 

M N T i-1 

P AVE (n »+ (,)) “ (n » l 7 ( J < t ,(s)ds) X N N 

1=1 ^ T i» T i-l ) 

where Xj denotes the characteristic function for the interval I 


defined 


(see [1] 


for details). 
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Corollary 4.5: There exist a costant K _> 0 (Independent of N ) 

such that for all z € Z 


IIP^ ve (P N [A N *] -1 P N z) - [A*] *zl < (*) |z|. 


Proof : A direct application of Theorem 4.4 yields the identities 


N 

\ (a 2k+l + a 2k ) = . 0 < 2k < N-l 


and 


N 


T <a Li + 4> ■ .1" - *J). 2 < a < «. 

i=2k 


On the other hand, since (?,4>(*)) - [A*] 1 (n ,4» (* ) ) - Lemma 4.3 implies that 


NT 0 
C = a Q = A (n + / <|>(s)ds) 

-r 


and 


s s 

>Ks) = (a^ + / A T (e)d9)5 - / <j> (e )de 

-r -r 


Therefore, if s^ = (tj + Tj_ 1 >/2, then 


I <*»1 + a " k >' 2 - 

N 

S 2 k+1 „ T 
= I / < B 2k+l (s) " 1)(A (sH ~ ^ s >) ds 

T 2k+1 


N 

T 2k 

+ / (B 2k+1 (s)(A (s)C - 4> ( s ) ) d s | 

S 2k+1 
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N N 

T 2k T 2k 

< (r//TT N) {(/ | A(s) | 2 ds) 1/2 |C | + (/ |*(s)| 2 ds) 1/2 }, 


N 
2k+l 


N 

x 2k+l 


and similarly, 




N 
2k- 1 


N 
2k- 1 


< (r//12 N) {( / | A(s) 1 2 ds) 1/2 |C | + ( / |<|>(s)| 2 ds) 1/2 } 


N 

2k 


N 

2k 


It follows that 


P N i N a 
*AVE 1 a 


N 


(?> l ^(s^)x 
j = l 


j ' ,v r N N 

lT J * T J-1 > 


)l 


< (r// 12N) (IAI, |§| + H<t»a T ). 

L 2 L 2 


Since f£»' ,2 (-r,0; If 1 ) and 


$(s) = a 1 (s) 5 - <j>(s) = a t (s)a T (n + / <j>(e)de) - <t>(s). 


-r 


there exists a constant K such that 


i* - Z *<s")x „ „ l, < (K/s)i(n,*(-))i z . 
J-l ] [tJ.tJ.,) L 2 - 


This estimate combined with the previous inequality establishes the proof. 


We turn now to providing a proof that the approximating adjoint semi- 


groups constructed above do not converge strongly to the adjoint semigroup 
generated by A*. 
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5. NON-STRONG CONVERGENCE 

Let A, A*, A^, and A^* be defined by (2.4) - (2.5), (2.6) - (2.7), 
(2.16) and (4.1), respectively. The corresponding semigroups will be denoted 
by S(t), S*(t), S N (t), and S^*(t). Recall that for z€Z (see [2]) 

(5.1) <A N z,z> - <AP N z,P N z> _< (ullP N zll 2 < toil zH 2 

where w = (1 + 2 1 Aq | + | A ^ ) 2 + 211 All ^ )/2 and for t _> 0 

(5.2) US N (t)ll < e ut , IIS N *(t)ll < e tot . 

The following result is a special case of Theorem 4.2 of Chapter 3 in [6]. 
Theorem 5.1: The following are equivalent: 

(a) For every z € Z and X € t with Re X > w 

N* -1 - * -l 

(XI - A ) z (XI - A ) z, as N + «. 

(b) For every z € Z and t ^ 0 

S N *(t)z + S*(t)z, as N + +°° , 

the convergence being uniform in t on bounded intervals. 


We shall also need the following technical lemmas 
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Lemma 5.2: Suppose the condition (b) of Theorem 5.1 holds. 


X £ p (A ) 


, i, r>N/-v T A N*v-l n N„ 
and HP (XI - A ) P II 


Is uniformly bounded in N. Then for every z € Z 


P N (XI - A'VVz * 


(XI - A*) 1 Z. 


Proof : From Theorem 5.1, for Xq > w and z€ Z 


P N (X 0 ! - 


.N*x-1 N ^ 
A ) P z 


<V - 


* ~1 
A ) z, 


as N -► «>. 


Note that for z € Z 


P N (XI - A N *) -1 P N Z 


- pS <V - 


A'VVz 


- (x 0 - X)P N (XI - a"*) 1 P N <X q I - A N *) _1 P N z, 


and similarly 


(XI - A*) *(I + (X - X Q )(X 0 I - A*) *) = (X Q I - A*) -1 . 


Hence, if w = z + (X - Xq)(XqI - A ) *z, then 


„N,. T ,N*,*1 N , 1T A *v— 1 

P (XI - A ) Pw-(XI-A) w = 


(X Q - X)P N (XI - A N *) 1 P N ((X q I - A N *) _1 P N z - (X Q I - A*) _1 z) 


+ P N (X Q I - A N *) -1 P N z - (X Q I - A*) _1 z. 
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This implies that for all z € Z 


P N (XI - A 8 Ww * (XI - A*)' 1 .. 


^ —1 

But since I + (X - Aq)(XqI - A ) is onto, the above statement holds for 
all w € Z. 

We turn now to a special case where ^ + al, a 6 r, 

A(s) = 0, -r £ s _< 0, A = Ag + and denote by A*, A N * and 

S (t), s (t) the corresponding infinitesimal generators and semigroups. 

Lena 5.3: If statement (b) of Theorem 5.1 holds for S N *(t), then it 

holds for < §^*(t). 


Proof: Note that 


* ~N* 

A - A + E and A 


, N* 


+ p n ep n 


where E: Z -*• Z is the bounded linear operator defined by 


E(n »<!>(•)) = (an, -A (*)n). 


V?JL 

It follows from (5.2) that II s ( t ) II _< e 0>c where Z - w + |a| + t! All 
Consequently, if X > u> » then 


u>t 


(XI - A N *) -1 = (XI - A N *) -1 + (XI - A N *) -1 P N EP N (XI - #*) 1 


(5.3) 
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and 

(5. A) (XI - X *) -1 = (XI - A*) -1 + (XI - A*) -1 E(XI - 'X*) -1 . 

Theorem 5,1 implies that for all z€Z 

(XI - A N *) _1 2 -► (XI - A*) -1 z 

Jc 

and since the ranks of E and E are finite it follows that 

N N 

IIP EP - Ell 0. 

It now follows from (5.3) - (5.4) that 

~ N * ^ — 1 , -1 

(XI - A ) z -► (XI - A ) z 

for all z€Z and this completes the proof. 

By Lemma 5.3, without loss of generality, one can assume that A(») = 0 

and A = Aq + Aj is invertible in what follows. We will show that there 

exists an element z£Z such that P^(A^*)”^P^z does not converge to 
*\-l 

(A ) z. First we consider the case when Aj is not the indentity. From 
Lemma 4.3, if (5,^(»)) = [A*] - *(ri,0) where 0 * n € I?, then 

5 = A T p and ip ( s ) = A^5 , -r _< s jC 0. 
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Applying Theorem 4*4 we obtain 


where 


pV'V'pVo) - iV 


N 

a 2k = S » 0 _< 2k _< N 


a 2k+l " + 2A 1 5 . 0 < 2k + 1 < N. 


For illustration, we have the following picture for the case N = 4 and n * 1 



where the solid line stands for (A*)”*(ti ,0) and the dashed line for i^a^ 1 . 
Since a^ k - A^£ » A^5 - a^k+l = 5 “ = e e Uf 1 (independent of N), it is 

easy to show that 


(5.5) IP N (A N *)" 1 P N (t 1 ,0) - (A*) -1 (n ,0)n| - 2N ||2. s e| 2 ds - I| e| 2 * 0. 

0 


Next we consider the case AJ = I. Let <|>(s) = x * 0 (constant vector in 
R n ). Then, from Lemma 4.3 (£ , ip(»)) - (A*) - ^ (0 ,<f> (• )) is given by 


£ = rA ^x and i|>(s) = £ - (s + r)x, -r < x < 0. 


(5.6) 
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And, from Theorem 4.4, P N (A N *) -1 P N (0 ,<J> (• ) ) = i N a N , where 


N r 

a o ? 


(5.7) 




l (£)x = 5 + ■£ rx, 0 < 2k < N 

jodd 
l<j<2k-l 


N 

a 2k+l 5 


2N-1 

N 


rx 


+ 2 l (J)* 

jeven 

2<j<2k 


,2N-2k-lv n . y „ 

= C - ( ^ )rx, 0 £ 2k < N. 


Since (a^ + a 2k+i^ 2 ^ T 2k + T 2k+1^ 2 ^ and ^ a 2k-l + a 2k^ 2 
^(( T 2k 1 + T 2k^^* ^°^ ows ^ rom (5*6) and (5*7) that 


r 

IN 2 

(5.8) IIP N (A N *)- 1 P N (0^(O) “ (A*) _1 (0,<(>(» ))\\y- 2N / |2Nsxj 2 ds - L- |x| 2 * 0. 


Now we may state the main theorem. 


Theorem 5.4: There exists an element z € Z and t > 0 such that 

S N *(t)z does not converge to S*(t)z. 

Proof : If for every z € Z and t > 0, S^*(t)z converges to S*(t)z, 

then it follows from Theorem 4.4 and Lemma 5.2 that for every z€Z 
pN(A^*) - lpNz (A*) - *z, where by Lemma 5.3 one can assume that A(») = 0 
and A = Aq + A^ is invertible. This contradicts the facts (5.5) and 


(5.8) 
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